Do Observations Favour Galileon Over Quintessence? 
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We study the Galileon scalar field model arising as a decoupling limit of the Dvali-Gababdaze- 
Porrati (DGP) construction for the late time acceleration of the universe. The model has one 
extra Galileon correction term over and above the standard kinetic and potential energy terms 
for a canonical quintessence field. We aim to study whether the current observational data can 
distinguish between this Galileon field and the quintessence field. Our study shows the remarkable 
result that for potentials like linear, square or exponential, the data prefers the Galileon model 
over quintessence field. It confirms that the observable universe demands the inclusion of higher 
derivative Galileon corrections in the standard quintessence scalar field models. 



Cosmological observations [lj [H indicate that our 
present universe is going through an accelerated expand- 
ing phase. The standard lore is that an unknown form 
of energy, called the dark energy 0], is responsible for 
driving the universe into such a late time accelerating 
phase. However, the nature of dark energy is still un- 
explained and remains a challange for particle physicists 
and cosmologists alike. So far, the simplest candidate for 
dark energy has been the cosmological constant which 
although is allowed by all observational data, yet it is 
plagued by acute problems like fine tuning and cosmic 
coincidence @. However, current observational data can 
also accommodate a time varying vacuum energy. Infact, 
quintessence @ ( a scalar field which slow rolls at present 
energy scale) was proposed as a candidate for dark en- 
ergy to provide a dynamical solution to the cosmological 
constant problem. 

Recently a large scale modification has been proposed 
which can explain the late time acceleration at the cos- 
mological scale. This involves an effective scalar field 7r 
dubbed as "Galileon" Q as its lagrangian respects the 
shift symmetry in the Minkowski background: n — > n + c 
and O^tt — > d^ir+b^ where c and are constants. Such a 
field can arise in the decoupling limit of the DGP model 
The lagrangian for such field can usually contain 
three terms: one linear in tt, one contains the usual ki- 
netic term for a canonical scalar field and the third one 
contains term like (V7r) 2 D7r. This third term in particu- 
lar is related to the decoupling limit of DGP model 
One can add another two terms involving higher deriva- 
tives in such a way that the final equation of motion for 
the ir field is still second order [Toj . This set up is theo- 
retically appealing due to the absence of negative energy 
instability as well as due to the absence of curvature sin- 
gularity. The model gives rise to late time acceleration 
of the universe 11| and at the same time is consistent 
in astrophysical context through Vainshtcin mechanism 

& 

Although theoretically Galileon model is more exciting 
(l3j than the standard quintessence scenario, the question 



is whether current observations can distinguish these two 
scenario. We aim to study precisely this question. Our 
goal is to see whether observations prefer the Galileon 
field over the standard quintessence scenario. 

We consider the action for the Galileon field in the 
lowest nontrivial order keeping upto the third order term 
in the lagrangian ( This arises in the decoupling limit of 
DGP model). We also keep a general potential term V(ir) 
in the action. 



S = 



rM n 2 , l 



where 



(1) 

8ttG is the reduced Planck mass, a is a 
dimcnsionlcss constant; for a = this action reduces to 
that of a standard quintessence field. V(ir) is the poten- 
tial for the 7T field. For V(tt) = ciir, it is the usual third 
order action for the Galileon field. S m is the action for 
the matter field. M is a constant of mass dimesnion one; 
by a redefinition of the parameter a, we can fix AI = M p \. 

Variation of the action ([1]) with respect to the met- 
ric tensor and assuming a flat Friedman-Robertson- 
Walkcr (FRW) spacetime with scale factor a(t), we get 
the Einstein's equations: 



3Af pl tf 2 = Pm + — ( 1 - 6^#7r ) + V(n) , 



(2) 



Ml (2H + 3H 2 ) = -- 1 + 2-3-tt + V(n) (3) 
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Varying the action ([1]) w.r.t the field it, we get the 
equation of motion for the field tt as 

(4) 

These above equations are supplemented by the matter 
conservation equation given by: 

p m + 3Hp m = 0. (5) 

Let us introduce the following dimensionless quantities 
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FIG. 1: Evolution of the equation of state uj^ for the ir field as function of redshift. From top to bottom, A,: = 1,0.5,0.1. 
f^mo = 0.3, e; = 15. The left one is for linear Potential and the right one is for exponential potential 



Then we have the autonomous system of equations: 

7T y/V 



(6) 



V6HM pl ' V3HM pl 
-6^^, A = -M pl ^ (7) 



3a; 3 (2 + 5e + e 2 ) - 3a; (2 - e + y 2 {2 + 3e)) + 2VEy 2 X - VEx 2 y 2 eX 

4 + 4e + x 2 e 2 

y (12 (-1 + y 2 ) (1 + e) - 6a; 2 (2 + 4e + e 2 ) + V&x 3 e 2 X + 2y/&x (2 + (2 + y 2 ) e) A) 

8 + 8e + 2a; 2 e 2 

e (-3a; (-3 + y 2 ) (2 + e) + 3a; 3 (2 + 3e + e 2 ) - 2^y 2 X - v^VeA) 



a; (4 + 4e + x 2 e 2 ) 



X' = V6a;A 2 (l-r), 



(8) 
(9) 
(10) 
(11) 



with T = i-^sp. and fi m = 1 - a; 2 (l + e) - y 2 . For e = 
we recover the autonomous system for the standard 
quintessence scenario flil ]. 

The equation of state for the n field is given by: 

-12y 2 (l + e) + 3x 2 (4 + 8e + e 2 ) - 2y/6xy 2 eX 
^ = 3(4 + 4e + a; 2 e 2 )(y 2 + a; 2 (l + e)) ' 

(12) 

In our system of equations (9)-(12), we have four vari- 
ables e.g x, y, e and A. The variable e sets the relative 
strength of the Galileon correction term over the standard 
kinetic energy term for a quintessence field. To solve this 
system, we need to specify the initial values for these four 
variables. We set our initial condition at the decoupling 
time z = 1000 and evolve the system from then onwards 
till the present epoch z = 0. 



Initially the it field is nearly frozen due to large hubble 
damping. This sets Xi to be small initially which keeps 
equation of state very close to —1. We vary Xi between 
0.01 and 10~ 7 and check that initial value of lu^ does not 
vary much and remains extremely close to — 1 . The initial 
value for e, e,, sets the initial strength of the Galileon 
correction over the the standard quintessence term. If 
this value is zero initially, it remains zero throughout the 
history of the universe and the evolution is same as the 
standard quintessence scenario. Hence the value of Ci 
determines the deviation from the quintessence scenario 
and the effect of the Galileon correction. So in our study 
we keep this parameter e, as a model parameter. The 
parameter A is the slope of the potential and its initial 
value determines the slope at which the field starts rolling 
initially. For smaller values of A,;, the ir field behaves 
very close to the cosmological constant (C.C). This is 
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FIG. 2: Evolution of the equation of state ui^ for the n field as function of redshift. Dashed, Dotted and Dash-Dotted lines 
represent ei = 0, 20, 40 whereas solid line represent the quintessence field. Potential is chosen to be linear one. The left one is 
for tlmo = 0.3 and the right one is for f2 m o = 0.24 



irrespective of the form of potential or whether the it field 
behaves as quintessence or Galilcon. As one increases A^ , 
the 7r field starts behaving differently from C.C. This is 
shown in Figure 1 where we plot the equation of state 
for the 7r field, ui w , for two potentials (Exponential and 
Linear) for different values of A,; . We plot this for Galileon 
field assuming e, = 15. 

Our aim is to distinguish between the quintessence and 
the Galileon field and hence we concentrate in the region 
where they deviate substantially from C.C. For this we set 
A, = 1 in our subsequent calculations which ensures that 
the 7r field deviates from the C.C behaviour at present. 

With the above choices for Xi, Ci and Ai, the initial value 
of the variable y, yi, is related to the present day matter 
density f2 m o. The parameter V controls the shape of the 
potential. In our study, we consider four different forms 
for the potential, e.g linear, squared, exponential and in- 
verse squared. These are the most well studied forms of 
potential for the scalar field model of dark energy. 

In Figure 2, we show the behaviour of the equation of 
state as a function of redshift for different values of 
ej. We choose two values for £l m0 e.g 0.24 and 0.3. We 
show it for the linear potential, but the overall behaviour 
remains the same for other potentials. From this figure, 
it is apparent that for = initially, the it behaves ex- 
actly same as quintessence. As one increases the ej, the 
two models starts deviating from each other. This devi- 
ation is slightly higher for smaller value of O m o- Hence 
constraining ej by observational data is crucial to distin- 
guish between quintessence and Galilcon field. 

In the rest of the paper, we concentrate on this issue. 
For this, we consider various observational data currently 
available. 

We consider the Supernovae Type la observation which 
is one of the direct probes for late time acceleration. We 
have utilized the Union2 compilation of the dataset which 
comprises of 557 datapoints |15j. It measures the appar- 



ent brightness of the Supernovae as observed by us which 
is related to the luminosity distance cIl{z) defined as 



d L (z) = (l + z) 



dz' 
o W) 



(13) 



With this we construct the distance modulus which 
is experimentally measured: 



/' 



M = 5 log 



Mpc 



25. 



(14) 



Here m and M are the apparent and absolute magnitudes 
of the Supernovae which are logarithmic measure of flux 
and luminosity respectively. 

Another observational probe that has been widely used 
in recent times to constrain dark energy models is related 
to the data from the Baryon Acoustic Oscillations mea- 
surements fpo| by the large scale galaxy survey. In this 
case, one needs to calculate the parameter D v which is 
related to the angular diameter distance as follows 



D„ 



zbao 

H(ZBAO) 



dz 
H(z) 



1/3 



(15) 



For BAO measurements we calculate the ratio 

j- ) t '| z _Q ■ This ratio is a relatively model indepen- 
dent quantity and has a measured value 1.736 ± 0.065. 

Next we use new determinations of the cosmic expan- 
sion history from red-envelope galaxies. Stern ct al. [TtJ 
have obtained a high-quality spectra with the Kcck-LRIS 
spectrograph of red-envelope galaxies in 24 galaxy clus- 
ters in the redshift range 0.2 < z < 1.0. They com- 
plemented these Keck spectra with high-quality, publicly 
available archival spectra from the SPICES and VVDS 
surveys. With this, they presented 12 measurements of 
the Hubble parameter H(z) at different redshift. The 
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measurement at z = was from HST Key project [l8j.At 
this point we define the normalised hubble parameter as 
h(z) = H ^ z ' and utilise it to derive the value of new h(z) 
and its deviation. 

Using all these observational data, our aim is to con- 
strain the parameter e; . We do not make any assumption 
on range of €4 which may arise due to the underlying 
theory. We simply want to see what is allowed by obser- 
vational data. 

The result is shown in Figure 3. Here we show the con- 
fidence contours in the q — S^mO parameter space for 
various potentials. We also show the bound on Cl m o 
as obtained by WMAP 7 observation which is given by 
n ro0 = 0.2669 ± 0.0288. The most striking result is 
that for both linear potential and square potentials, the 
quintessence case (e; = 0) is ruled out by more than 
95% ( 2a) confidence level. For the exponential poten- 
tial, the quintessence case is ruled out at 68% (la) confi- 
dence level. Only for the inverse squared potential both 
the quintessence and the Galileon behaviours are allowed 
both at la and 2a confidence level. We should stress that 
the linear potential is the most simple and well studied 
case for Galileon field. Hence for this simple potential as 
well as for other potentials, quintessence field is clearly 
distinguishable from the Galileon field by current obser- 
vational data. This result also confirms the fact that the 
observational data actually prefers the Galileon correc- 
tion term in the scalar field lagrangian. 



To conclude, the inclusion of Galileon correction term 
in the standard canonical scalar field lagrangian is well 
motivated and elegant. Usually one can include five terms 
in a scalar field lagrangian which respect Galileon sym- 
metry atleast in the flat case and also result second or- 
der equation of motion avoiding the appearances of in- 
stability. Our aim is to see whether such inclusion of 
higher derivative terms is also observationally interesting. 
In other words whether cosmological observations prefer 
such terms over the standard canonical kinetic term. We 
stress that our goal is not to compare Galileon model with 
ACDM which is consistent with all the current oberva- 
tions. Hence we consider the parameter region where the 
7r field behaves differently from C.C. 

We restricts ourselves to the simple case of Galileon 
model which contains three terms in the lagrangian and 
which arises in the decoupling limit of the DGP model. 
By fitting this model with observational data, we show 
that data indeed prefers this Galileon correction over the 
standard canonical kinetic term for most of the poten- 
tials including the most simple linear potential. This, in 
our knowledge, is the first result confirming that current 
observations demand the inclusion of Galileon correction 
in the scalar field lagrangian. This makes the Galileon 
model even more interesting. 
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